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ABSTRACT 

We present the first spectroscopic analysis of the faint M31 satellite galaxies, And XI 
and And XIII, as well as a reanalysis of existing spectroscopic data for two further 
faint companions, And IX (correcting for an error in earlier geometric modelling) and 
And XII. By combining data obtained using the DEep Imaging Multi-Object Spectro- 
graph (DEIMOS) mounted on the Keck II telescope with deep photometry from the 
Suprime-Cam instrument on Subaru, we have identified the most probable members 
for each of the satellites based on their radial velocities (precise to several km s" 1 down 
to i ~ 22). Using both the photometric and spectroscopic data, we have also calculated 
global properties for the dwarfs, such as systemic velocities, metallicities and half-light 
radii. We find each dwarf to be very metal poor ([Fe/H]~ —2 both photometrically 
and spectroscopically, from their stacked spectrum), and as such, they continue to 
follow the luminosity-metallicity relationship established with brighter dwarfs. We are 
unable to resolve dispersion for And XI due to small sample size and low S:N, but we 
set a one sigma upper limit of a v <4.5 kms -1 . For And IX, And XII and And XIII 
we resolve velocity dispersions of 0^=4. 5^34, 2. 6*2.® an d i^jj^J kms" 1 , though we 
note that the dispersion for And XIII is based on just three stars. We derive masses 
within the half light radii for these galaxies of 6.2±jj;? x 10 6 M Q , 2.4^1 x 10 6 M Q and 
LI^qy x 1O 7 M0 respectively. We discuss each satellite in the context of the Mateo 
relations for dwarf spheroidal gala xies, and in referenc e to the Universal halo profiles 
established for Milky Way dwarfs ( Walker et alll2009D Both And IX and XII fall be- 
low the universal halo profiles of IWalker et al.l (|2009T ) . indicating that they are less 
massive than would be expected for objects of thei r half-light radius. When combined 
with the findings of iMcConnachie fc Irwin! ()2006aD . which reveal that the M31 satel- 
lites are twice as extended (in terms of both half-light and tidal radii) as their Milky 
Way counterparts, these results suggest that the satellite population of the Andromeda 
system could inhabit halos that with regard to their central densities - are significantly 
different from those of the Milky Way. 

Key words: cosmology: dark matter - galaxies: abundances - galaxies: Local Group 
- galaxies: dwarfs - galaxies: individual - galaxies: kinematics 
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In recent years, the low luminosity population of dwarf 
spheroidal galaxies (dSphs) has drastically increased in num- 
ber, with ~30 faint satellites discovered both within the 
halo of our Milky Wa y (MW) and in that of our "sis- 



ter" g alaxy, M31 (e .g 



Irwin et al. 



2008 



20051; iMartin et all l200rj: 



Zucker et all |2004|; IWillman et all 

20071; 



.Zucker et all 2006b a 
McConnachie et al.ll2008l ; iBelok 



et al.l 



^ 

2008; Mar tin et all 120091 ). Whilst in principle these minus 



cule stellar associations should be some of the simplest ob- 
jects of the galactic family - which spans roughly seven or- 
ders of magnitude in mass - to understand, the analysis of 
faint dSphs (My > -8) within the MW has dissolved the 
notion that they form in a simple, homogeneous fashion. De- 
tailed studies (both photometric and spectroscopic) of such 
objects reveal several unusual properties for the population. 
For example, they highlight a significant departure from 
the well-estab lished mass-luminosity relationship for dSphs 
i|Mateolll998l ). A number of these newly discovered objects 
are also quite controversial in their nature, possessing prop- 
erties that are common to both dSphs and globular clus- 
ters and are therefore difficult to classify as either, such as 
the unusual objects, Willman 1 and Segue 1 (IWillman et al.l 
l2005l ; iBelokurov et al.ll2007l ; iNiederste-Ostholt et al.ll2009l ). 
In addition, comparisons of the MW and M31 dSph sys- 
tems have revealed significant differences between the pop- 
ulations, where the classical M31 dwarfs (My < —8) are 
observed to be at least twice as extended as their Milky 
Way counterparts, in terms of both their half-light and tidal 
radii (iMcConnachie fc Irw in 2006a). In addition, the spatial 
distribution of the M31 satellites around their host is more 
extended than the Galactic satellite distribution, with ap- 
proximately half of the Galactic satellites found within 100 
kpc of the Galaxy, compared to roughly 200 kpc fo r M31 
|McConnachie fc Irwinl [2006bl ; iKoch fc Grebelll2006l ). Such 
differences indicate that environment plays a significant role 
in the evolution of dSphs, the full extent of which is not 
known. 

Another significant point of interest concerning dwarf 
galaxies is the well known discrepancy between the ob- 
served number of these satellites, and the number pre- 
dicted by ACDM models in numer ical simulations (e.g. 
iMoore et allll999l ; iKlypin et al.lll999l '). which still amounts 
to one to two orders of magnitude. Survey completeness is 
likely to contribute t o at least part of this deficit (see e.g. 
Willman et all 120041 ; ISimqn fc Gehal 120071; llbllerud et al.1 



2008; iKoposov et al 



120081 ; IMcConnachie et al.ll2009l) 



the Galactic case, observational efforts are hampered by ob- 
scuration from the disk and the bulge. It is thought that 
future all-sky surveys may find a wealth of ultra-faint satel- 
lites that would reduce the gap between observation and 
theory. Other explanations assume that many of these ha- 
los would remain dark to the present day, owing to sup- 
pression of s tar formation as the result of a photo-ionizing 
background (|Somerv illc 2002). These halos would therefore 
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represent a significant proportion of the satellite halos found 
in simulations, reconciling these with observations. 

In an attempt to better understand this unusual popula- 
tion of galaxies, and to unite theory with observation, spec- 
troscopic surveys of faint dwarf galaxies (particularly within 
the MW via the Sloan Digital Sky Survey (SDSS)) have been 
carried out. Using the velocity information from the spec- 
tra, one can measure the systemic velocities and intrinsic 
velocity dispersions for these objects, whi ch can be used to 
estimate the mass of these systems (e.g. Illingworthl Il976l; 



20091 ; IWolf etail 



iRichstone fc Tremainelll986l ; IWalker et al 
2009). However, methods such as these make various as- 
sumptions about the anisotropy and virial equilibrium of 
the systems that are not necessarily correct. Nonetheless, 
such methods have been shown to be goo d indicators of 
the instantaneous mass of the system ( e.g. lOh et all 1 19951 ; 
iPiatek fc Prvorlll995l ; iMufioz et al.ll2008l ) , and can therefore 
provide an insight into the mass distribution of the satel- 
lites. Historically, the results of such studies indicated the 
existence of a mass limit of ~ 10 7 Mq under which no dwarf 
galaxies are found, supporting the theory that these satel- 
lites would indeed only inhabit massive dark matter halos. 
The existence of such a m ass lim i t wou ld confirm the em- 
pirical relation defined by iMateol l|l998h for satellites that 
are brighter than Draco (Mv = —8.8) and Ursa Minor 
(M v = -8.9) in the Local Group: M/L = 2.5 + (L/10 7 L Q ), 
where M/L is the mass-to-lig ht ratio of the galax ies and 
L is their luminosity. Work by IStrigari et all (2008D . inves- 
tigating the mass contained within the central 300 pc of 
dSphs, reinforced this relation, demonstrating that, despite 
the studied dSphs spanning five orders of magnitude in lu- 
minosity, they appear to show a lower dynamical mass limit 
of ~ 10 7 Mq . However, such a finding is in stark contrast 
with the masses derived for some of the recently discovered 
ultra-faint dSphs, e.g. Segue 1 and Leo V, whose masses of 
- 10 5 M Q (|Geha et all 120091 ; IWalker et all I2009D fall well 
below t his bound. Th i s disc repancy has been addressed in 
work by I Walker et~"ai1 (|2009l ), where they evaluate the mass 
of dSphs within the half-light radius (ri/2). Their analy- 
sis shows that these ultra faint dSphs may simply be more 
deeply embedded within dark matter halos that are similar 
to those inhabite d by their larger and bri ghter counterparts 
(also discussed in iPenarrubia et aflbOOSl ). which could sug- 
gest the existence o f a "universal" dar k matter halo. Similar 
results are shown in i Wolf et alj l|2009h . where they conclude 
that all the Milky Way dSphs are consistent with inhabit- 
ing dark matter halos with a mass Mhaio ~ 3 x 1O 9 M0. 
The recently discovered faint dSphs described above pop- 
ulate the faint end of the satellite luminosity function, and 
their spectroscopic analysis would allow us to see if these ob- 
jects are indeed highly dark matter dominated, and whether 
they conform to such a mass limit. 

Until recently, only a handful of the M31 dSphs have 
been analysed spectroscopically (And II, And I X, And XII, 
And XIV, And XV And X VI a n d And X JCote et al l 
19991: Ichapman et al.l l200fj|. 120071; iMaiewski et all 120071 



Letarte et al.1 120091 ; iKalirai et al.l 120091 ). with only one of 
these (And XII by Chapman et al. 2007) probing the fainter 
en d of the luminosity function (Mv >-8). A recent paper 
bv lKalirai et al.l (|2010l ) has extended the number of bright 
dSphs surveyed, presenting new data for And I, II and III, 
but the kinematic properties of M31's lowest luminosity 
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dSphs are still unknown. Spectroscopic observations of these 
objects would allow us to see how they compare to both their 
brighter M31 brethren, and their MW cousins of a compa- 
rable luminosity. 

As a step towards this end, we have initiated a kine- 
matic survey of th e trio of faint dwarfs discovered by 
iMartin et all l|200fj ). And XI, And XII and And XIII using 
the DEep Imaging Multi-Object Spectrograph (DEIMOS) 
on Keck II. We combine these kinematics with deep photom- 
etry from Suprime-Cam on the Subaru telescope. The anal- 
ysis of our observations is performed here for these objects, 
including a re-treatmen t of And XII, whose o rbital prop- 
erties were modelled in IChapman et all (j2007T h but whose 
detailed properties were not assessed. We also include a re- 
examination of the fa int satellite, And IX, u sing the same 
spectroscopic data as IChapman et all (|2005l ). as there was 
an error in the geometric modelling of the slit mask used to 
calculate the velocities, which led to a misclassification of 
member stars; this has since been corrected for, and we also 
present additional Subaru photometry for the satellite. The 
outline of the paper is as follows: § 2 presents the observing 
strategy, the observations and how they were reduced ; § 3 is 
dedicated to the analysis of And XI, And XII And XIII and 
And IX; we discuss our findings in § 4 and provide conclu- 
sions in § 5. 



2 OBSERVATIONS 

2.1 Photometric Observations 

Each of the satellites were observed photometrically using 
the Suprime-Cam instrument, mounted on the primary fo- 
cus of the 8-metre Subaru telescope, located at the summit 
of Mauna Kea. This instrument is ideal for probing the spa- 
tial extent of these faint dwarfs as it can reach the Horizon- 
tal Branch (HB) of populations at distances greater than 
500 kpc in a reasonable exposure time. Also, it has a wide 
field of view (34' x 27') which is more than sufficient to cover 
the entire region of the dwarfs at the distance of M31. 

Observations of And XI, And XII and And XIII were 
made on the night of September 22, 2006, in photometric 
conditions with an average seeing of 0.6" and 0.5"in the g 
and i bands respectively. For these satellites we exposed 
for a total of 2100s in g-band (split into 5 x 420s dithered 
subexposures) and 4800s in the i-band (split into 20 x 240s 
subexposures) in an attempt to reach below the HB to an 
equal depth in both bands. The data were processed using 
the CASU general purpose pipeline for p rocessing wide-field 
optical CCD data jlrwin fc Lewis! 120011 ). The images were 
debiased and trimmed before being flat-fielded and gain- 
corrected to a common internal system using master flats 
constructed from twilight sky observations. Catalogues were 
generated for every science image and used to refine the as- 
trometric alignment of the images. The images were then 
grouped for individual objects and passbands and stacked to 
form the final science images based on the updated astrom- 
etry. A catalogue was then generated for each final stacked 
science image, objects morphologically classified as stellar, 
non-stellar or noise-like, and the g,i band catalogue data 
merged. As a final step, for And XI, XII, XIII the g,i Subaru 
data was calibrated with respect to our MegaCam system 



g,i using all stellar objects in common with signal-to-noise 
> 10. And IX was observed in service mode during the pe- 
riod 9-10 October 2004 with typically subarcsecond seeing, 
for a total of 600s in each of g, r and i bands (split into 5 x 
120s dithered subexposures for each filter). The data were 
reduced in a similar manner to those for And XI, XII and 
XIII, although here, the photometry from the CASU Subaru 
pipeline was matched to the SDSS M31 scan data. 

The data for all satellites are extinc tion corrected and 
dereddened using the lSchleeel et al Il998l dust maps. We find 
E(B-V)=0.08 for AndXI, XIII and IX, and E(B-V)=0.12 for 
And XII. 



2.2 Spectroscopic Observation 

For each of the dSphs in question, candidate members were 
selected using colour magnitude diagrams (CMDs) from 
the Canada France Hawaii Telescope (CFHT) as shown 
in Fig. [T] where we present the MegaCam colours (with- 
out correcting for extinction and reddening) for the stars 
within these objects, transformed i nto Landolt V and I us- 
ing a two stage pro c ess d etailed in llbata et alj (|2007l ) and 
iMcConnachie et al.l |2004l ). Colour magnitude boxes were 
used to select RGB stars that lay within 3-5' regions sur- 
rounding each of the satellites, which is equivalent to sev- 
eral half-light radii of these objects. DEIMOS is a multi-slit, 
double-beam spectrograph mounted on the Nasmyth focus 
of the Keck II telescope on Mauna Kea, ideal for taking 
optical spectra of faint-objects such as these. For AndXI, 
And XII, and And XIII multi-object observations were taken 
from 21-24 September 2006 with DEIMOS in photometric 
conditions and with excellent 0.6" seeing. Our chosen instru- 
mental setting covered the observed wavelength range from 

o 

5600-9800 A with exposure time of 3x20 minute integra- 
tions. The Red Giant Branch (RGB) stars in these three new 
dwarfs are faint (i = 21.5-23.0), so to obtain a useful S:N in 
a reasonable exposure time, all the spectra employ the low 

o 

resolution (~ 3.5 A) 600 line/mm grating. AndIX was ob- 
served on the night of 13 September, 2004 under photomet- 
ric conditions with seeing from 0.5" to 0.8". We employed 
the 1200 line/mm grating, covering the wavelength range 

o 

6400-9000 A, also with exposure time of 3x20 minute inte- 
grations. Data red uction was perform ed using the DEIMOS- 
DEEP2 pipeline (|Faber et al'l \2Q0$ \, which included debi- 
assing, flat-fielding, extracting, wavelength-calibrating and 
sky subtracting the spectra. In table [T] we detail the posi- 
tions targeted by each mask, the number of target stars, and 
the total number of stars for which accurate velocities were 
derived (discussed further below). 

The kinematic data f or An dIX were originally pre- 
sented in IChapman et all (|2005l '). However the geometric 
model used for the slit mask was incorrect, causing a shift 
in the measured velocities of the observed stars of order 
5kms _1 , which has been corrected in this work. This system- 
atic error caused a misclassification of member stars within 
the sample, which in turn affected the measurements of sys- 
temic velocity, dispersion and metallicity for the satellite in 
Chapman et al. (2005). The velocities for all observed stars 
has since been re-measured allowing us to correctly derive 
the properties of And IX in this work. 
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Figure 1. CFHT-MegaCam color-magnitude diagrams (left panel of each Figure) , where the n a tive MegaCam g an d i bands have been 
converted into V and I magnitudes using transformations detailed inlMcConnachic et al. 

I l l2004) ;l Ibata et al. I 11200711 . and radial velocity 
uncertainties (right panel of each Figure) of the observed stars in the four objects surveyed in the paper. The top left panel corresponds 
to And XI, the top right to And XII and the bottom left and right correspond to And XIII and And IX respectively. The small dots 
represent CFHT objects in the region, filled circles represent target stars selected in the satellites CMD Red Giant Branch features. 
Circled big points represent stars that were initially selected as field stars but have the radial velocity and the CMD position of dwarf 
stars and thus were used to derive the dwarf parameters. 



Table 1. Details of the Keck/DEIMOS observations of AndXI, XII, XIII and IX 



Date observed 


Mask a./2ooo 


Mask 5j2000 


Object 


Grating 


Field P.A. 


Exp. time (s) 


No. targets 


No. reduced 


23/09/2006 


00:46:21 


+33:48:21 


AndXI 


600 line/mm 


0° 


3x1200 


35 


32 


21/09/2006 


00:47:27 


+34:22:30 


And XII 


600 line/mm 


0° 


3x1200 


50 


48 


23/09/2006 


00:47:27 


+34:22:30 


And XII 


600 line/mm 


0° 


3x1200 


50 


48 


23/09/2006 


00:51:51 


+33:00:16 


And XIII 


600 line/mm 


0° 


3x1200 


49 


48 


13/09/2004 


00:52:51 


+43:11:48 


And IX 


1200 line/mm 


0° 


3x1200 


648 


277 



2.2.1 Velocity calculations and error analysis 

Our instrumental setup covered the wavelength range of 
the Calcium II (Ca II) triplet, a strong absorption feature 

o 

centered at ~ 8550 A. To determine velocities for our tar- 
geted stars, we cross-correlated each spectrum with a model 
template of the feature at the rest-frame wavelength posi- 
tions of the triplet lines (the technique used is similar to 
that discussed by Wilkinson et al. 2004). Briefly, we per- 
form cross-correlations for each star with our template and 



compute the corrected heliocentric velocity for each star. 
We also compute the Tonry-Davis R-value (TDR) for each 
cross-correlation, which is an indication of the reliability 
of the resulting velocity. We consider any velocity cross- 
correlation with a TDR value>3.5 to be reliable. In Ta- 
bles 2 - 5, we note the TDR values for each of our likely 
member stars (discussed in greater detail in §3). We also 
estimated the velocities for each star using a Doppler shift 
method, fitting Gaussians to each of the triplet lines indi- 
vidually and computing the mean velocity shift from these. 
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Figure 2. Telluric velocity corrections for all stars in And XII as 
a function of I— band magnitude (top) and S:N (bottom). The 
median telluric correction for this mask is 2.8 kms -1 , and for 

o 

bright (I < 21), high S:N (S:N> 12 A" 1 ) targets, this correction 
remains approximately constant. However for fainter stars with 
noisier spectra, this correction shows increased scatter about the 
median value. 



Figure 4. The total error in velocity of all stars observed within 
And XI, XII and XIII (top panel) and And IX (lower panel) as a 
function of S:N. The total error is calculated by combining the 
Monte-Carlo error (ctmc) discussed in §2.2.1, and the system- 
atic error, e (6.2 kms" 1 for And XI, XII and XIII, 2.2 kms" 1 for 
And IX), in quadrature. It can be seen that we can measure the 
velocities of our sample to accuracies of typically less than 15 

kms" 1 down to S:N>3 A" 1 . 




-2 2 

( 7 B =(v I -v 2 )/((75 + tr|+2E 2 ) 1 / 2 (km/s) 

Figure 3. A histogram showing the normalized error distribu- 
tion for repeat measurements of the same stars in our And XII 
data set. The normalized error, crjv incorporates the velocity dif- 
ferences between the repeat measurements (vi and vj), and the 
Monte-Carlo uncertainties calculated for each observation (ci and 
(T2). In order to reproduce a unit Gaussian distribution for our 
uncertainties, we also include an additional error term, e, which 
accounts for any systematic uncertainties we have not included. 
For the And XII data (observed with the 6001 grating), we find 
e=6.2 kms -1 . For And IX, which was observed with the 12001 
grating, we set 6=2.2 kms -1 , which was the value derived by 
ISimon &; Gehal (120071 ) from repeat observations of stars in their 
Keck/DEIMOS survey, which used a similar observational setup 
to that of And IX. 



The results of both calculations were generally found to 
be consistent (within ~ 2 — 5kms -1 ) with one another. 
In a few instances, the Gaussian technique revealed cases 
where the cross-correlation technique had fixed onto a noise- 
feature, rather than the true Ca II triplet, and these cases, 
we have taken the velocities derived from the Gaussian 
technique. However, as it is generally accepted that sim- 
ple Doppler-shift velocitie s are less reliable in low S:N data 
(e.g. lAllende Prietoll2007h . we refer to velocities computed 
by the cross-correlation method throughout the text, unless 
otherwise stated. 

For the Ca II triplet lines, the spectral features of inter- 
est overlap in location with strong night sky OH recombina- 
tion lines for some of the dSph members. While we are able 

o 

to detect continuum at S:N>3 A -1 for even the faintest 
(i — 23) stars observed, the velocity errors are subject to 
larger systematic errors when they lie close to these OH 
lines. This effect is seen most prominently in the case of 
And XII, where the Call8498 line coincides with a strong sky 
feature. This effect is seen in the night to night variations 
in the velocities of the And XII stars, for which we obtained 
repeat observations. A systematic shift in the velocities be- 
tween the two nights of 30 kms -1 is observed for some of the 
stars, while other stars in the mask at more favourable ve- 
locities have night to night errors <5 kms -1 . For those stars 
where this effect is very prominent, we mask out the affected 
line and repeat the cross-correlation on the two unobscured 
lines. 

There are a number of systematic errors that can be 
introduced into velocity calculations from the instrument 
setup, such as the effect of slight mask rotation, field mis- 
alignment and improperly cut masks, all of which can cause 
the star to be mis-centred within the slit, resulting in an 
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incorrect velocity measurement that can range from a few 
to tens of kms -1 . To correct for this error, we make use of 
telluric absorption features present within our spectra, par- 
ticularly the strong Fraunhofer A-band telluric absorption 

o 

feature located between 7600 and 7630 A. These features are 
the result of absorption of stellar light by various gases (O2, 
H2O, CO2, etc.) within the Earth's atmosphere. These lines 
are therefore not Doppler shifted, and should be present at 
the same wavelength in each sp ectru m. Following the tech- 
niques of ISimon fc Gehal (120071 ) and lSohn et~al] (|2007h . we 
construct a template telluric spectrum which can be cross- 
correlated with each of our science spectra to determine the 
effective velocity shift caused by these systematic errors. For 
the fields observed in this work, we found that the average 
velocity shift for each mask derived by this method ranged 
from +2.8 kms -1 to +8.1 kms -1 . To correct our observed 
velocities for this, we add the velocity shift calculated for 
each star to our cross-correlation velocities determined from 
the Ca II triplet. In Fig[2]we show the telluric corrections 
derived for And XII (average correction = 2.8 kms -1 ) as a 
function of /—band magnitude and S:N. For brighter stars 

o 

(I < 21) and stars with higher S:N (S:N> 15 A" 1 ), the ve- 
locity corrections are approximately constant, however for 
fainter stars (stars with lower S:N) the magnitude of this 
telluric correction shows significant scatter about the me- 
dian value. It is therefore unclear whether this method can 
accurately compute the telluric shift in fainter objects. 

We estimate the errors on o u r velo cities by following 
the pr ocedures of ISimon fc Gehal (|2007l ) and iKalirai et al.l 
|20ld ). First, we make an estimate of our velocity uncer- 
tainties for each observed star using a Monte Carlo method, 
whereby noise is randomly added to each pixel in the spec- 
trum, assuming a Poisson distribution for the noise. We then 
recalculate the velocity and telluric offset of each star, us- 
ing the same cross-correlation method detailed above. This 
procedure is repeated 1000 times, and then the error is cal- 
culated to be the square root of the variance of the result- 
ing mean velocity. Secondly, we make use of our And XII 
dataset, where we were able to obtain repeat observations 
for our entire mask, giving us 47 stars with two velocity mea- 
surements, v\ and «2- We then attempt to measure our veloc- 
ity errors directly by comparing these velocities and defin- 
ing a normalised error, gm which incorporates our Monte 
Carlo error estimates for each measurement pair, 01 and 02, 
plus a random error, e, which represents any uncertainties 
we have not accounted for. e is defined as the additional 
random error required in order to reproduce a unit Gaus- 
sian distribution with our data (shown in Fig. [3]). For the 
And XII stars, which were observed with the 600 line/mm 
grating, we find e — 6.2 kms" 1 . While this is roughly double 
that found in other work s (e.g. 2.2 kms -1 inlSimon fc Gehal 
120071 and ~ 3kms -1 in IKalirai et alJl201Ch . we note that 
they use the 1200 line/mm grating, which has a resolution 

o 

of ~ 1.3 A compared to the 600 line/mm gratings resolution 

o 

of ~3.5 A. Therefore, for Andr omeda's XI, XII and X III, we 
set e = 6.2 kms -1 , and we use ISimon fc Gehal (|2007l ) value 
of e = 2.2 kms -1 for And IX. We then combine this system- 
atic error in quadrature with our calculated errors from the 
Monte Carlo technique (cmc), giving us our final error esti- 
mate. We display the errors for all observed stars in the four 



dSphs as a function of S:N in Fig. [4] The dSphs observed 
using the 600 line/mm grating (And XI, XII and XIII) are 
shown in the top panel, and And IX (observed with the 1200 
line/mm grating) is shown in the lower panel. It can be seen 
in both cases that we can measure velocities to accuracies of 
~ 5— 15kms -1 down to a S:N of 3 or better, while the higher 
signal-to-noise stars have velocity errors that are consistent 
with the 2.2 and 6 kms -1 systematic uncertainties for the 
1200 and 600 line/mm gratings respectively. We find a me- 
dian error of atot = 3.6 kms -1 for stars in the AndIX mask, 
and atot = 6.6 kms" 1 for the stars in the And XI, XII and 
XIII masks. 



3 SATELLITE PROPERTIES 
3.1 And XI 

And XI is the most luminous of the three dSphs reported in 
iMartin et al. I (|2006l l. It was found as an over-density of stars 
in the CFHT/MegaCam survey that follows a well-defined 
RGB. The satellite lies at a projected distance of 106 kpc 
from M31 (assuming a distance of 785 kpc for M31). Below 
we present the structural and kinematic properties of this 
faint galaxy. 

3.1.1 Photometric properties 

First, we discuss the updated structural parameters derived 
for the And XI dwarf using the data from Subaru/ Suprime- 
Cam. The extinction corrected and deredenned CMD for 
this data is shown in Fig. [5{a), containing only stars 
within a radius of 1.5' of the centre of And XI, defined 
as the mean of the positions of the stars within 3', giving 
(a,5) = (00 h 46 m 21 s , +33° 48' 22"). The Subaru data probe 
much further down t he luminosity functi on than the original 
CFHT data used bv lMartin et all (120061 ) and shown in Fig. 
1, and is complete down to the level of the HB at g~25. To 
assess the half-light and tidal radii of And XI, we construct a 
background corrected radial surface density profile of And XI 
using circular annuli; this can be seen in Figure[5jc). A back- 
ground level of 8.1+0.4 stars per square arc minute was de- 
termined from a series of comparison regions in the area sur- 
rounding the dwarf. The half-light radius is then calculated 
by using a chi-squared fitting routine to apply King (with 
rt = 2!23 ± 0.2), Plummer and exponential profiles to the 
radial density data. For And XI this yields half-light radii 
(r 1/2 ) of 0!64±0.15, 0!72±0.06 and 0!71±0.03 respectively, 
with the exponential profile providing the best fit (with a re- 
duced chi-squared o f 2.6) to the data. A lternatively, apply- 
ing the algorithm of IMartin et al.l (|2008l ) produces a similar 
result of 0!68±0.04 and converges on structural parameters 
without a strong ellipticity, 77 = 0.08 (N. Martin, private 
communication), justifying the use of circular annuli in this 

case. 

The CFHT data analysed bv lMartin et all (|200eJ) were 
not deep enough to measure the distance to And XI, how- 
ever our deeper, more complete Subaru observations allow 
us to place constraints on the distance to the satellite by 
combining the photometric data with our spectroscopic data 
(which we discuss in more depth below). One can define the 
apparent magnitude of the Tip of the Red Giant Branch 
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Figure 5. Top left: Here we display the de- reddened Subaru CMD for And XI and highlight the stars for which DEIMOS spectroscopy 
is available with mage nta squares. The ou tlier at -432.3 kms -1 is highlighted in green. A Dartmouth isochrone with [oj/Fe]=+0.4,[Fe/H= 
-2.0 and age 13 Gyrs (Dottcr ct al. 2008) is overlaid at a distance modulus of 24.5 (equivalent to 760 kpc) with [Fe/H]=-2.0. Top right: 
The luminosity function for And XI in the g-band is displayed as a histogram, and a clear peak at 25.1 for the HB is observed compared 
with the normalised luminosity function of a comparison region in the surrounding area of the dwarf (blue histogram). Bottom panel: 
The radial density plot for And XI with King (solid line), Plummer (long dashed line) and exponential (short dashed line) overlaid, with 
corresponding r 1/2 = 0^64+0.05, 0^72+0.06 and 0(71+0.03. 



Table 2. Properties of candidate member stars in And XI, centered at a=00:46:21, 5=+33:48:22. Outlier at v = — 432.3kms 1 is also 
included as the last entry. 



star a 


star <5 


vel (kms 1 ) 


[Fe/H] spec 


[Fe/H] phot 


S:N ( A" 1 ) 


TDR 


/-mag 


V-mag 


00:46:22.65 


+33:46:21.6 


-418.8+6.2 


-2.42 


-2.0 


16.7 


15.37 


20.5 


21.9 


00:46:19.10 


+33:48:4.1 


-417.2+6.3 


-1.56 


-1.5 


13.5 


17.2 


20.8 


22.3 


00:46:15.14 


+33:48:26.7 


-419.8+6.2 


-0.61 


-1.9 


9.3 


8.4 


21.3 


22.5 


00:46:21.56 


+33:48:22.4 


-421.9+6.4 


-1.69 


-1.4 


7.8 


7.3 


21.5 


22.9 


0:46:16.96 


+33:48:3.8 


-432.3+10.8 




-1.8 


6.0 


3.4 


21.8 


22.9 
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Figure 6. Left panel: The velocities of observed stars in the And XI field are shown as a histogram, with the 5 candidate And XI 
member stars (seen as a kinematic grouping at approximately -420 kms^ 1 ) highlighted as a heavy histogram. To differentiate AndXI 
stars from the field, we additionally plot the velocities against their radius from the And XI center, and in the lower panel, we display their 
photo metric metallicities as a function of radial velocity, which were derived by interpolating between Dartmouth isochrones ijDotter et al] 
2008) with [«/Fe]=+0.4 and age=13 Gyrs. Right panel: The normalised spectra for each of the 4 plausible member stars of AndXI 
and the outlier (bottom spectrum) are shown, with the combined spectrum of the 4 plausible members plotted in the top panel. The 
spectroscopic [Fe/H] for each candidate member is shown also, and we derive [Fe/H]=-2.0±0.2 from the composite spectrum. 



(TRGB) for the dwarf by assuming that our brightest spec- 
troscopically confirmed member represents the tip, allow- 
ing us to calculate the distance to the satellite. Using the 
brightest DEIMOS star (which has an apparent magnitude, 
Io,Ve 9 a=20.393), and assuming t he absolute magnitud e of 
the TRGB is M; 8 , Ve9 a = -4.04 jBellazzini et al.ll2004l ), we 
find the maximum distance to be T) max = 770 kpc. To deter- 
mine the minimum distance, we must assess what the max- 
imum possible offset could be from this brightest RGB star 
to the true TRGB. T o do this, we follow the procedure of 
IChapman et al.l l|2007t ) , analy sing the well populated RGB o f 
the dwarf spheroidal, And II (|McConnachie &: Irwinll2006al ) . 
With 1000 random samples of 4 of the brightest stars from 
the top magnitude of the And II RGB, we constrain a prob- 
ability distribution for the true TRGB, corresponding to a 
most likely offset of 0.5 mag. Using this value, we find the 
range of distances to AndXI is 610-770 kpc. We note that 
this calculated range is based on assumption that the bright- 
est DEIMOS star is found within 0.5 mags of the TRGB. 
While this seems reasonable, the distance range we calcu- 
late would vary if this assumption were changed. Therefore, 
as a sanity check for these values, we can also look for the 
peak in the luminosity function due to the HB of And XI in 
the luminosity function, shown in Fig. [5jb). An inspection 
of this reveals a clear peak at g-magnitude of ~ 25.1 (red 
dashed line). Given that the absolute magnitu de of the HB 
in the g-band is M ws =0.6 (|lrwin et al.ll2007l ) - though we 
note that this as sumed value is se nsitive to age and metal- 
licity effects, see IChen et al.ll2009l for a discussion - we can 
estimate a distance modulus for the satellite of ~ 24.5 ( 830 
kpc), which sits outside of our calculated range for the satel- 
lite. This result demonstrates the difficulty in determining 



distances for faint satellites such as AndXI. A companion 
paper, analysing HST data for AndXI, XII and XIII (Mar- 
tin et al. in prep), will discuss the distances to these galaxies 
in more detail. For now, we use an intermediate distance of 
760^150 kpc (where the upper and lower error bounds cor- 
respond to the maximum and minimum distances estimated 
within this work) for the remainder of our analysis. 

In order to assess a mean photometric metallicity for 
the satellite, we us e the Dartmouth isochrones models of 
iDotter et al.l ([2008), at an alpha element abundance of 
[a/Fe]=+0.4 and an age of 13 Gyrs. We use [a/Fe1=+0.4 as 
it ha s been shown in numerous studies (e.g. iMarcolini et al.l 
2008), that the alpha abundance for metal-poor stars in 
dSph populations such as those studied herein, tend to scat- 
ter between [a/Fe]=+0.2 and +0.5. We overlay isochrones 
at distances of 610 kpc, 830 kpc and 760 kpc (correspond- 
ing to the upper, lower and HB distances calculated above), 
and note the best observed fit in each case. This gives a 
metallicity range of [Fe/H]=-1.5 to -2.0 for the dwarf. In 
Fig. [S] we overlay the isochrone corresponding to the best 
observed fit at 760 kpc ([Fe/H]=-2.0). We note that using 
a less alpha-enhanced isochrone model (with [a/Fe]=-|-0.2), 
gives a slightly more metal rich value, with a best observed 
fit at the HB distance of [Fe/H]=-1.9. 

We use the Subaru data to recalculate the luminosity 
of the satellite by summing the flux of the likely member 
stars lying on the RGB and HB in the CMD of the satel- 
lite. By summing the flux of stars populating the CMD of 
the well studied globular cluster M92_J within its half light 



We use M92 as a comparison, rather than another M31 dSph 
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Figure 7. Likelihood distributions of member stars in And XI, showing the most likely values of radial velocity (left panel) and velocity 
dispersion (right panel) for our 4 confirmed members. The dashed lines indicate the 1, 2, 3xcr uncertainties. The most-likely velocity for 
And XI is v r =-419.6 kms -1 and the dispersion is not quite resolved. We use the value of 4.6 kms -1 from the la confidence level as an 
upper limit for the dispersion of And XI. 



radius (r 1/2 =19.6" iFerraro et afl|200d ). using archival HST 
data, we calculate the total flux contained in this radius 
from the TRGB down to the HB. As our own data start 
come incomplete at roughly the HB level, we calculate a 
range for the contained flux by choosing two different cuts 
for the HB of M92. We find that 58-68% of the total lumi- 
nosity of the satellite is contained within this region, and 
we calculate a value for the total luminosity of And XI of 
4.9^q'2 x 1O 4 L0. Using the revised distance estimate above 
and this luminosity, we update the derived parameters of 
And XI, with the half-light radius, T1/2, becoming 145^29 
pc using the half-light radius from the exponential profile 
(previously estimated as 115 pc), and M v = —6. 91° ^, in 
agreement with the findings of iMartin et al.l (|2006l ). 

3.1.2 Spectroscopic properties 

Here we present the spectroscopic results from And XI; the 
data for each spectroscopic member star can be found in Ta- 
ble [2] Of the stars that were identified as potential And XI 
members from the CFHT CMD and targeted with DEIMOS, 

as its CMD is well populated past the HB level, allowing us to 
better estimate the fraction of luminosity contained within this 
region. In addition, the satellites studied in this paper all appear 
to be relatively old, making the comparison with a globular cluster 
reasonable. 



four appear to be kinematically associated with the satellite. 
As shown in the top left panel of Fig. [6l where we out- 
line the spectroscopic results for And XI, one can clearly 
identify a kinematic peak lying at a heliocentric velocity 
of ~ — 420kms _1 . These four stars have a mean velocity 
of v r = — 419.4kms _1 and an unusually low dispersion of 
o~ v — 2.4kms _1 (which we note is uncorrected for instru- 
mental errors). We also initially identify one other star as a 
potential member of AndXI with a velocity of -432.3 kms -1 . 
While it is clearly an outlier in velocity, it is important to as- 
sess its potential membership to the satellite. The inclusion 
of this star as a member decreases the uncorrected values of 
v r to -422 kms -1 and a v to 4.0 kms -1 , however the velocity 
error on this fifth star is greater than the other 4 candidate 
members (10.8 kms -1 vs. average v err of 6.2 kms -1 ). This is 
partially due to the fact that the CaIIs662 line is completely 
buried in the noise of the spectrum, as shown in the bottom 
right panel of Fig. [6] As a result, when cross-correlating the 
spectrum with the Call template we achieve a peak in the 
correlation function of less than 0.1 and a TDR value of only 
3.4, implying that our velocity measurement for this star is 
unreliable. For this reason, we define this star as a tentative 
member of AndXI, but we exclude it from the following 
kinematic analysis. 

It is important for us to determine whether any of our 
remaining stars are contaminants in the form of either fore- 
ground G dwarf stars or M31 halo stars. Because AndXI 
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Figure 8. As Fig.[5]but for And XII. A Dartmouth isochrone is overplotted with [Fe/H]=-1.9 and [a/Fe]=+0.4 at distance modulus of 
24.6 (830 kpc). In plot (c) Plummer and exponential profiles are overlaid, with respective half-light radii of 1'.2 ± 0.2 and 1' 1 ± 0.2, with 
the exponential profile providing the best fit to the data. . A King profile (with rj/2 fixed at 1.2' is overplotted also. 



Table 3. Properties of candidate member stars in And XII, centered at a=00:47:27, <5=+34:22:29. 



star a 


star 5 


vel (kms 1 ) 


[Fe/H] spe c 


[Fe/H] pfcot 


S:N ( A- 1 ) 


TDR 


/-mag 


V-mag 


00:47:31.04 


+34:24:12.1 


-566.9+6.6 


-2.4 


-2.0 


8.9 


5.0 


21.1 


22.4 


00:47:24.69 


+34:22:23.9 


-567.8+6.3 


-2.3 


-1.7 


10.1 


10.7 


21.2 


22.5 


00:47:27.76 


+34:22:6.2 


-549.8+6.2 


-1.7 


-1.7 


11.4 


11.0 


21.2 


22.6 


00:47:28.63 


+34 22 43.1 


-561.2+6.3 


-2.1 


-1.8 


9.7 


8.9 


21.5 


22.8 


00:47:31.34 


+34 22 57.6 


-549.3+7.0 


-1.9 


-1.6 


7.3 


6.75 


21.7 


22.9 


00:47:26.65 


+34 23 22.8 


-547.4+6.5 


-2.5 


-1.7 


6.2 


5.97 


22.0 


23.3 


00:47:27.18 


+34 23 53.5 


-557.8+6.7 


-2.0 


-1.7 


5.4 


5.3 


22.0 


23.2 


00:47:30.60 


+34 24 20.3 


-557.4+7.3 


-1.9 


-1.5 


4.8 


3.7 


22.2 


23.3 



has such a highly negative velocity (~ — 420kms _1 ), con- 
tamination from MW G dwarfs is unlikely, as these are pre- 
dominantly found at velocities of gr eater than -160 kms -1 
l|Robin et al.ll2004l : flbata et all 120051 ). so the chances of us 
finding these stars at the velocity of And XI are slim. Never- 



theless, we can make use of the Na I doublet, a pressure sen- 

o 

sitive absorption feature observed at ~ 8100 A. The strength 
of this feature is strongly correlated with the surface gravity 
of stars, therefore it is prominent in G dwarfs, and much 
weaker in RGB stars. As our spectra cover the region of this 
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Figure 9. As Fig. [6] but for And XII. In the left panel, And XII can be clearly identified as a kinematic grouping of 8 stars at approximately 
-560 kms , all located within 2' of the dwarf centre. The individual stars are all found to be metal poor, both photometrically (lower 
left panel) and spectroscopically (right panel). Analysis of the composite spectra of the confirmed members yields [Fe/H]=-2.1±0.2 
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Figure 10. Likelihood distributions for And XII. Dashed lines represent the conventional 1,2 and 3a confidence levels. The mean velocity 
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for And XII is calculated to be -558.4 kms 1 , and we (marginally) resolve a velocity dispersion of 2.6^j g kms 1 , though we note that 



this is consistent with zero within its lcr uncertainties. 
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Figure 11. As Fig. [5] but for And XIII. Dartmouth isochrones jDotter et al.ll2008h is overlaid with [Fe/H]=-2.0 and [a/Fe]=+0.4 at 
the HB distance modulus of 24.8 (910 kpc as seen in the luminosity function of AndXIII, plot (b)) In plot (c), King (with it = 2.8'), 
Plummer and exponential profiles are overlaid, with respective half-light radii of 0'78 ± 0.08, 0'81 ± 0.09 and 0(76 ± 0.1, with the King 
profile providing the best fit to the data. 



doublet, we can measure the EWs of the lines. Any star for 

o 

which the EW of the Na I doublet is greater than 1.8 A we 
classify as a foreground contaminant. As there is some over- 
lap in the EWs of these two populations, this will inevitably 
lead to the occasional misclassification, but it should en- 
sure that the majority of foreground stars are removed. For 
the four stars we classify as members of And XI we mea- 

o 

sure the EW of the Na I doublet to be less than 1.5 A in 
all cases. Contaminating M31 halo stars are more difficult 
to cleanly define. The M31 halo has a systemic velocity of 
~ — 300kms~ 1 and an average dispersion of ~ 90kms _1 , 
|Chapman et alj|2006l ; iKalirai et afll2006h so it is possible 
that stray halo stars in the tails of the distribution could 
be found at a velocity similar to that of And XI. However, 
given the position of And XI, found in a low-density region 
of M31's halo, finding such a star within a few half-light 
radii of the object is unlikely. To limit any such contami- 
nation, we impose a radial cut on our members, requiring 



them to sit within the tidal radius of the dSph, 2.13'. We 
plot the distance of all observed stars from the centre of 
And XI as a function of their velocity in Fig. [6] It can be 
seen that all our candidate members lie within this radius. 
As a final check for contaminants, we plot the photomet- 
ric metallicity of each star as a function of velocity, derived 
by interpolating be tween the same Dartmouth isochrones 
|Dotter et alj 120081) as used in § 3.1.1. Halo stars are ex- 
pected to have metallicities of [Fe/H]~ —1.4, and foreground 
dwarfs analysed with these isochrones would likely appear to 
have metal-rich metallicities of [Fe/H]~ —0.6. Our 4 mem- 
ber stars and the outlier are all found to be similarly metal 
poor, with a median [Fe/H]=-1.9. Combining these factors, 
we find no reason to reject any of our four plausible members 
as contaminants. 

The spectra for each of the candidates (including the 
outlier) are displayed in Fig. [6] along with a composite spec- 
trum for the four most plausible members (weighted by their 
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Figure 12. As Fig. [6] but for And XIII. In the left panel, And XIII can be identified as a kinematic grouping of 3 stars at approximately 
-200 kms" 1 , all of which are located within 2' of the dwarf centre. The individual stars arc all found to be metal poor, both photometrically 
(lower left panel) and spectroscopically (right panel). Analysis of the composite spectra of the confirmed members yields [Fe/H]=-1.9±0.2 




Figure 13. Likelihood distributions for And XIII. Dashed lines represent the conventional 1,2 and 3cr confidence levels. The mean velocity 
for And XIII is determined to be -195.0 kms , and the velocity dispersion is resolved as a v = 9.7±|;§ kms" 1 . 
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signal-to-noise), which is shown in the top right panel, high- 
lighted in blue. As each of the spectra has a high enough S:N, 
we can determine their metallicity individually and also cal- 
culate an average value from their sum, f rom the Calcium II 
triplet (as described in llbata et al.|[2005l ). giving a resulting 
measure of [Fe/H] for the satellite which is independent of 
photometric calculations. This relation between the EWs of 
the Ca II triplet and the resulting metallicity is dependent 
on the apparent V-band magnitude of the observed star. By 
summing together stars of different magnitudes we introduce 
an error into our calculations as we must assume an aver- 
age value for the V-band magnitude of our sample. As the 
stars analysed within this work are all of a similar luminos- 
ity, this error is typically small (less than 0.1 dex), and we 
factor this into our error calculations. The average value de- 
termined from the combined spectrum is [Fe/H]= — 2.0±0.2. 
T his is consider a bly m ore metal poor than that calculated 
bv lMartin etall l|2006h of [Fe/H] =-1.3 from the CFHT pho - 
tometry. This is likely due to the use of lGirardi et all (|2004T ) 
isochrones in their study, which produce more metal-rich val- 
ues in the MegaCam filter sets than other isochrone models, 
such as the Dotter isochrones used in this paper. This effect 
has been noted in the glo bular clusters NGC 2419 and M92 
dGirardi et al .2004,2008), where the Girardi et al. isochrone 
models in the g-i colours produce more metal-rich values for 
these objects than other published values. 

Following the ma ximum likelihood algorithm outlined 
m iMartin et all (|2007l ). we calculate the mean radial veloc- 
ity, v r , and the intrinsic velocity dispersion, a, for And XI 
by sampling a coarse grid in (v r ,a) space and determining 
the parameter values that maximise the likelihood function 
(ML), defined as: 



ML(v r 



exp 



O"tot 



(1) 



velocity measured for the i 
uncertainty and otot 



star, 



the radial 
r ,i the corresponding 
. t . In this way, we are 
able to separate the intrinsic dispersion of And XI from the 
dispersion introduced by our measurement uncertainties. We 
display the one dimensional likelihood distributions for v r 
and a v in Fig. [7] where the dashed lines represent the con- 
ventional 1, 2 and 3a (68%, 95% and 99.7%) uncertainties 
on the values. We recover a systemic velocity for And XI 
of u r =-419.6l3 g kms -1 (where the errors indicate the la 
uncertainties). Owing to our low sample size and S:N, the 
dispersion is unresolved, however the la value suggests a 
maximum value of 4.6 kms -1 from the data, which we over- 
lay as a Gaussian curve in the top panel of Fig. [6^,. We note 
that these values are only marginally affected if we include 
the tentative member at -432.3. kms -1 , with v r increasing 
to -422.5 kms -1 , and a v is still unresolved at kms -1 , with 
a maximum set by the la uncertainty of ~ 5kms -1 . 

These results suggest that the velocity dispersion of 
And XI is unusually low, comparable to the curious MW 
objects Willman 1 and Leo V, which have measured 
disper sions of 4.2+ 2 ' 1 kms -1 and 2.4+ 2,4 kms -1 respec- 



tively jMartin et al.ll2007l ; I Walker et alj|2009h . and the M31 



dSph, And X, wh ich has a dispersion of 3.9±1.2kms 
jKalirai et al.ll2009h . We caution that the calculation of dis- 
persions for cold systems such as And XI rely heavily on 



accurately knowing the errors attributed to each velocity 
measurement. Whilst we have attempted to investigate our 
errors as thoroughly as possible (as detailed in §2.2.1), if 
we have overestimated our errors, we will have underesti- 
mated our velocity dispersion (and vice- versa), and the true 
velocity dispersion could be higher than the dispersion we 
calculate here. 

Using the spectral data obtained above, we can make an 
estimate for the mass of And XI by assuming it is a spherical 
system in virial equilibrium. Throughout the literature, sev- 
eral different methods for calculating the mass of dispersion- 
supported galaxies - such as And XI- are employed which 
do not necessarily produce similar values when applied to 
the same data sets. This makes comparisons of these values 
for Local Group dSphs difficult, as one is not necessarily 
comparing like with like. In this work, we shall use three 
different methods to e s timat e the mass of our dSp hs; the 
iRichstone fc Tremaind (|l986h and llllingworthl (|l976l ) meth- 
ods, which are the two most common approaches, and w e 
also apply a new method, presented in lWalker et al.l ([2009) , 
which we discuss in more detail below. 

First, we appl y the core-fitting method of 
IRichstone fc Tremaind (Il986l ). This method assumes 
that the system is in virial equilibrium with a flat core, an 
isotropic velocity distribution and a constant mass-to-light 
ratio throughout. One can then derive the mass-to- light 
ratio (M/L) of And XI from: 



M/L = T) 
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2ttG Sorhb 



(2) 



where 5*o is the central surface brightness of the dwarf, r^b is 
the value of the radius where the surface brightness drops to 
half the peak surface brightness (approximately 2/3 x r 1( / 2 ) 
and (Jo its central velocity dispers ion. r\ is a model de- 
pende nt, dimensionless constant that lRichstone fc Tremainel 
( 1986) showed to be very nearly unity for a wide range of 
mass models, e.g. the isothermal sphere, where they found 
77=1.013. For And XI, we derive So = 0.455 ± 0.09 L pc -2 
when using fj,v = 27.35 mag/arcsec 2 , and taking ri/2 = 
145^29 P c we fi n d r hb = 97^19 pc. We further assume <jq — 
a v = 4. 7 kms -1 (taking the la limit from our maximum 
likelihood calculation) and use this to derive an effective up- 
per limit on the mass-to-light ratio for And XL This leads 
to (M/L) < 152Mq/Lq. Combining this with the updated 
luminosity of 4.91q2 x 10 4 L© calculated above, we deduce 
an upper limit on the mass of And XI of M < 7 .4 x 10 6 M . 

The prescription of llllingworthl (Il976l ). uses an 
r 1//4 law to describe the surface brightness profile of 
the system, rather than t he core-fitting technique of 
IRichstone fc Tremaind (l986). It again assumes virial equi- 
librium and an isotropic velocity distribution for the system. 
One can then derive the mass using: 



M 



(3) 



where 77=85OM pc -1 km -2 s 2 , we calculate M<2.6xl0 6 M 
and M/L=53, which is less than half the value derived above, 
highlighting the difficulty faced in comparing the masses of 

dSphs within the Local Group. 

Recently, papers bv lWalker etall (|2009h and lWolf et al.l 
( 2009) have shown that the mass within the half-light radius 
of dispersion-supported galaxies, such as dSphs, can be ac- 
curately estimated with only mild assumptions about the 
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spatial variation of the stellar velocity dispersion anisotropy 
(/3), giving it an advantage over other methods (such as those 
described above), which are more sen sitive to the adopte d 
value of p. Following the approach of I Walker et al.l (2009), 
we can derive the mass contained within r 1 / 2 (M-half) from: 

Mhalf = (J,r 1 / 2 a r v (4) 

where fi = 580Mopc -1 km -2 s 2 . For AndXI, this gives us an 
upper limit on M half for the dSph of 1.7xl0 6 M o . We do 
note that, as 2 of our 4 AndXI members sit outside of the 
half-light radius for the dwarf, we are assuming that there is 
no variation of a v with radius, which may not be the case. 
This could cause us to underestimate the mass of AndXI 
within this region. 

3.2 And XII 

And XII was the fa intest of the trio of dwarfs presented in 
iMartin et all £2006), with an absolute magnitude of — - 
6.4 . Its orbital properties were discussed in a previous paper 
by IChapman et ail (|2007l ). where it was shown that due to 
its large radial velocity, it was likely only experiencing the 
gravitational effects of M31 for the first time. In this sec- 
tion, we give a complete overview of the photometric and 
spectroscopic properties of And XII in the same format as 
for AndXI. We list the properties of our spectroscopically 
confirmed members in Table [3] 

3.2.1 Photometric properties 

First, we discuss the structural properties derived for 
And XII from the Subaru data. We determine the centre of 
AndXIIas (a,<5) = (00 h 47 m 27 s , +34° 22' 29") using the same 
technique as for AndXI, and we display the Subaru CMD 
for all stars within a 1.5' radius of the centre of And XII 
in Fig. [HI a) (pink squares represent the DEIMOS candidate 
members in this data set). The background-corrected radial 
profile can be seen in Figure (HJc), with a background level 
of 9.9±0.4 stars per square arc minute. Plummer and ex- 
ponential models were fitted and are shown over-plotted on 
the data. The Plummer and exponential fits yield half-light 
radii of l!2±0.2 and l!l±0.2 respectively, with the plummer 
profile showing the best fit to the data (reduced chi-squared 
of 3.7 vs. 3.9). The King profile would not converge on a 
physical fit for these data, yielding results where the tidal 
radius is over an order of magnitude smaller than the half- 
light radius. We therefore fit a King profile, with the half 
light radius fixed to the half-light radius given by the Plum- 
mer profile. The resulting best fit tidal radius is r t =3.l'±0.3, 
and we use this value as a cut for membership in §3.2.2. As 
can be seen in Fig. [8l the radial profile for And XII is not 
as clean as for AndXI, and this may be due to the pres- 
ence of variable extinction in the field of And XII, which 
appears (based on visible cirrus) to vary on angular scales 
of ~ 1', and is reflected by the larger errors in (of order 
~ 10 — 20% compared with ~ 5% for AndXI) and the inabil- 
ity to successfully fit a King profile to the data. This small 
scale ex tinction is also pres ent in a nearby globular cluster, 
MGC1 (|Mackev et al.ll2009h . and is found to vary on scales 
much smaller than the 6' x 6' resolution of the Schlegel et al. 
(1998) dust maps of this region. Conversely, these poor fits 



could be indicative of tidal disruption in the outer regions 
of And XII; however, given the large systemic velocity of - 
558.4 kms -1 (discussed in § 3.2.2), this is unlikely to have 
been as a result of an interaction between And XII and M31, 
as And XII is l ikely experiencing the potential of its host for 
the first time (|Chapman et al.ll2007l ). 

Adopting the same TRGB Monte Carlo procedure as 
with AndXI we find D maa; = l Mpc, and assuming a max- 
imum offset of 0.5 mag from the true TRGB, we compute 
a range of distances to And XII of 800-1000 kpc, consistent 
with the most likely distance f or And XII of 830 ±50 kpc 
found bv lChapman et all l|2007h . An inspection of the lumi- 
nosity function of And XII shown in the top right panel of 
Figure |8jb) reveals a clear peak in the HB at g-magnitudc 
just below ~ 25.3, implying a distance modulus for the satel- 
lite of 24.7 ( 870 kpc), which sits within our calculated 
range. We use t h e mos t-likely distance of 830^30° kpc from 
IChapman et al.l l|2007h for the remainder of our analysis. 
Again, we recalculate the luminosity of the satellite, find- 
ing a value of L = 3.i1q' 2 x 1O 4 L0, corresponding to an 
absolute magnitude of M„=-6.41q 5 and a half-light radius 
(using the value from the exponential model), T1/2, becomes 
289^49 P c (previously estimated as 125 pc). 

We analyse the CMD of And XII to ascertain a mean 
value of [Fe/H] for the dSph. We overlav lDotter et al.l (|2008T ) 
isochrones (with age and abundance as detailed in § 3.1.1) 
for the range of distances calculated for And XII above (800- 
1000 kpc), and conclude that the metallicity of And XII lies 
between [Fe/H] =-1.7 and [Fe/H]=-2.0. In Fig [8] we overlay 
the best observed fit for the most-likely distance (830 kpc), 
with [Fe/H]=-1.9. 

3.2.2 Spectroscopic properties 

We now discuss the kinematic results for And XII. As seen 
in Fig O 8 stars are found to be kinematically associated 
with And XII at a velocity of ~ —560 kms -1 . Unlike AndXI 
we have no ambiguity about which of these stars are mem- 
bers as they are tightly distributed in velocity, all sit within 
2'(< 2r!/ 2 ) of the dwarf centre, and all show negligible ab- 
sorption at the position of the Na I doublet. Their photo- 
metric metallicities (lower left panel of Fig. |5J) are also in 
good agreement with one another, ranging from [Fe/H]=- 
1.5 to -2.0, with median [Fe/H]=-1.7. We display the indi- 
vidual spectra for each candidate member in Fig. [§1 with 
the combined spectrum for the satellite shown in the top 
panel. Using the Call approach as described for AndXI, 
metallicity measurements are made for each of the candidate 
stars, spanning [Fe/H]=-1.5 to -2.5, with a median [Fe/H]=- 
2.0. To estimate an average metallicity for the dwarf, we 
analyse the composite spectrum, and deduce a metallicity 
of [Fe/H] = —2.1 ± 0.2. Eliminating the 3 stars with poor 
S:N (bottom 3 panels of Fig. |9j) yields a value of [Fe/H]=- 
2.0±0.25. Both of these values agree well with the photo- 
metric metallicities discussed above. 

Before correcting for measurement errors, the 8 stars 
imply a v = 4.9 kms -1 . The same maximum likelihood ap- 
proach as with AndXI yields a systemic velocity for the 
satellite of -558. 4^3'?, kms -1 , and we (marginally) resolve a 
velocity dispersion for the satellite of a v = 2.6I2 6 kms -1 , 
where the errors represent the la uncertainties. Again, this 
value is uncommonly low for a dSph. 
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Table 4. Properties of candidate member stars in And XIII, centered at a=00:51:51, <5=+33:00:16. 



star a 


star S 


vel (kms 1 ) 


[Fe/H] spec 


i^/H] phot 


S:N ( A -1 ) 


TDR 


/-mag 


V-mag 


00:51:53.00 


+33:00:18.2 


-202.6+6.2 


-1.9 


-1.7 


9.9 


16.2 


21.0 


22.4 


00:51:49.81 


+33:00:41.1 


-195.7+6.5 


-1.9 


-1.6 


3.3 


7.2 


22.1 


23.3 


00:51:42.34 


+32:58:05.3 


-184.7+6.6 


-1.8 


-1.9 


4.7 


5.9 


21.9 


23.1 



As before, we apply eqns. (2)- (4) to our And XII data in 
an attempt to constrai n the mass and dark matter c ontent 
of the galaxy. Using the lRichstone Sz Tremainel (1 19861 ) equa- 
tion for mass-to-light (eqn. (2)) and our updated parameters 
for r 1/2 of 289±% pc, S = 0.2l±° ° l L pc -2 and the disper- 
sion value from the maximum likelihood analysis of S.Slg'g 
kms -1 , we find M/L— 55^° 7 Mq/ Lq, giving a mas s esti- 
mate of l.Ttl.7 x 1O 6 M . Eqn. (3) (|lllingworthlll976r ) gives 
us M/L=54±^| 7 and a mass of \.7±\l x 10^ M B , in agree- 
ment with the iRichstone fc Tremainel (|l986h approach. In 
both cases, this value falls below the historical ~ 10 7 Mq 
mass threshold typically observed for brighter dwarfs, even 
when taking the large errors into account. 

Finally, applying the same formula for the mass within 
the half-light radi us of And XII as for And XI (eqn. (4) 
I Walker et al]|2009l ), we calculate M half = l.lt?;J x 10 6 M o . 
For And XII, six of our 8 member stars lie within the half- 
light radius of the dSph, making our measured a v a reason- 
able approximation for the dispersion within Ti/2- 

3.3 And XIII 

Discovered bv lMartin et al.1 ((2006), And XIII suffered some 
chance contamination with a background galaxy cluster, in- 
creasing the apparent contrast of this dwarf in the imagery. 
Nonetheless, the apparent brightness of And XIII is simi- 
lar to And XI after correction for this contamination. We 
present our observations of And XIII in the same format as 
the other M31 dwarfs, and list our spectroscopic members 
in Table H 

3.3.1 Photometric properties 

Using the same Subaru photometric data set as for 
And XI and And XII, we find the centre of And XIII as 
(a,5) = (00 h 51 m 50 s , +33° 00' 43") and in Fig. [Til; a) we show 
the Subaru CMD for all stars within a 1.5' radius of the 
centre of And XIII (DEIMOS candidate members are high- 
lighted with pink squares, while tentative members are high- 
lighted by light blue squares) . The background-corrected ra- 
dial profile is displayed here also (panel (c), with a back- 
ground level of 11.2+0.4 stars per square arc minute). King 
(with r t = 2.8 ± 0.2 ), Plummer and exponential models are 
shown overplotted on the data. These yield half-light radii 
of 0!78 ± 0.08, 0!81 ± 0.09 and 0!66 ±0.1 respectively, with 
the King profile showing the best fit to the data. Again, this 
value of r-i/2 is consistent with th e rj/g resulting from the 
application of the ML algorithm of iMartin et alJ |200&t ) (N. 
Martin, private communication). 

Adopting the same TRGB Monte Carlo procedure as 
with the previous 2 dwarfs, we find D mal =940 kpc for the 
satellite, and assuming a maximum offset of 0.5 mag from 



the TRGB, a range of distances of 750-940 kpc. We also 
find a distance modulus from the HB of And XII of 24.8 
(910 kpc), via inspection of the luminosity function, which is 
in good agreement with the results from the TRGB method. 
Using the HB distance, we calculate a total luminosity for 
the satellite of L = 4.lto'i, x 10 4 L Q . These revised dis- 
tance and luminosity values update the derived parame- 
ters of And XIII with M„ becoming -6.7jlg '\ (previously es- 
timated as -6.9), and the half-light radius, r^j, becomes 

pc (previously esti mated as 115 p c). 

As before, we overlay iDotter et al.l ((2008) isochrones 
(with age and abundance as detailed in § 3.1) for the range of 
distances calculated for And XII above (750-940 kpc), giving 
a metallicity range of [Fe/H]=-1.8 to -2.1 for AndXIII. In 
Fig [H] we overlay the best observed fit for the HB distance 
(910 kpc) which has [Fe/H]=-2.0. 

3.3.2 Spectroscopic properties 

Being the most distant of our 4 observed dSphs, AndXIII 
has unsurprisingly provided us with the biggest challenge 
in identifying probable member stars. When examining the 
individual spectra for the stars in the AndXIII mask (a sam- 
ple of which can be seen in Fig. I12[) we see that the CaIIg662 
line has experienced contamination from a skyline, which 
has significantly affected its profile. As a result, the EW of 
the third line, which should be ~ 3/4 the EW of the sec- 
ond line, is larger than that of the second in most cases. 
Such an effect will obviously change the measured veloci- 
ties and metallicities if we were to include it in our cross- 
correlations and metallicity calculations. For this reason, we 
derive velocities and metallicities for the satellite members 
based on the Call8498 and CaIIs542 lines only. The resulting 
kinematic data for And XIII is shown in Fig. 1121 We can see 
that three RGB stars lying on the Subaru CMD (Fig. [TTJ 
inhabit a similar region in velocity space, with an average 
v r ~-190 kms -1 . This implies that AndXIII is currently 
moving away from M31. These three stars have an uncor- 
rected dispersion of cr„=4.3 kms . The systemic velocity 
of AndXIII places it much closer to the regime of Galactic 
foreground contamination (vuel > —150 kms -1 ) than the 
previous two dwarfs. With this being the case, it is partic- 
ularly important for us to assess whether any of our spec- 
troscopic candidates could be contaminants. Measuring the 
EW of the Nal doublet in each star, we find that in each case 

o 

the EW is significantly less than the 1.8 A cut we defined 
in §3.1.2, making them unlikely foreground dwarf contami- 
nants. To assess the likelihood of any of the three stars being 
M31 halo stars that happen to be coincident with AndXIII 
we invoke a radial cut equal to the tidal radius (2.8'). All our 
potential member stars lie within this cut. Finally we assess 
the quality of our derived velocities, and we note that all 
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Figure 14. As Fig.[5]but for And IX. Isochrones are overlaid with [Fe/H]=-2.2 at the HB distance modulus of 24.5 (760 kpc). In plot (c), 
King (with rt = 6.3), Plummer and exponential profiles are overlaid, with respective half-light radii of 2'.7 ± 0.2, 2(6 ± 0.1 and 2(5 ± 0.1. 



three stars have TDR> 3.5, meaning that they satisfy our 
quality criteria. Therefore we use all three potential mem- 
bers in our subsequent analysis. Using the summed spectrum 
of these members we calculate a metallicity for And XIII 
of [Fe/H]=-2.0±0.3, which is a gain more metal po or than 
that of [Fe/H]=-1.4 reported in Martin et all (|2006l ), but in 
good agreement with our revised photometric metallicities 
detailed above. 

The maximum likelihood distribution for the veloc- 
ity and dispersion of And XIII is shown in Fig. 1131 
For our two probable members, the technique sug- 
gests v r =-195.0tg'4 kms -1 and a resolved dispersion of 



7+8.9 



kms . The corresponding errors are large (> 
50%) owing to our small sample size (3 stars) and the low 

S:N (9.9, 4.7 and 3.3 AT 1 ) of our objects. If we com- 



pare this value with dispersions for MW objects of a com- 
par able luminosity such as CVnl (tr„=7.9 kms~ , My=- 
7.1. iMartin et al.N2008l ; ISimon fc Geha|l2007l ) and Bootes I 
(<r„=6.5 kms~ 1 .Mv=-6.3. lMartin et al.ll2008l 120071) . we find 
that the dispersion of And XIII is consistent with its coun- 
terparts within its la errors. 

We combine these results with the structural proper- 
ties of And XIII derived from the Subaru data, to place a 
constraint the dark matter content of And XIII, using the 
same approach as for the previous dwarfs. Using eqn. (2) 
and a central surface brightness of So=0.276^q'q|Lq pc -2 
and the dispersion from our maximum likelihood analy- 
sis of a v — 9.7 kms -1 , we calculate a mass to light ra- 
tio of M/L = 82ll^f Mq/L q which implies a high dark 
matter dominance. The upper limit for the mass is there- 
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Figure 15. As Fig. [6] but for And IX. AndlXcan clearly be identified as a kinematic grouping at ^-200 kms -1 . Of the nine stars at 
this velocity, 8 are associated with And IX, with the 9th star lying well off the dwarf RGB (see Fig. 1141 at a distance greater than 10 
arcmins from the dwarf, and the two stars with low quality cross-correlations (TDR<3.5) highlighted in pink. The remaining 8 stars are 
all found to be metal poor both photometrically (lower left panel) and spectroscopically (right panel), however we note that 2 stars have 
unreliable velocity measurements. An analysis of the composite spectrum of the confirmed members yields [Fe/H]=-2.2±0.2. 
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Figure 16. Likelihood distributions for And IX. Dashed lines represent the conventional 1,2 and 3<r confidence levels. The mean velocity 
for And IX is calculated to be -207.7 kms -1 , and the velocity dispersion is resolved by the technique as 4.5 kms -1 . 
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fore 3.4±2.2 x 10 7 M Q . The formula of llllingworthl ||l976h 
(eqn. (3)), produces an upper limit on the mass of And XIII 
of I.6+10 x 1q7 M (corresponding to M/L=390 + 5 5 2 05 ), which 
is less than half the va l ue de rived using the method 
of iRichstone fc Tremaind jl986h . as was the case for 
And XI. Applying the formula for M^aif as before, we find 
Mhalf =1. 1±H x 1O 7 M . 



3.4 And IX 



First reported by the SDSS team in IZucker et all (|2004h . 
And IX has previously been studied both photometrically 
JZucker et al.ll2004l ; lHarbeck etldl2005l;lMcConnachie et all 
l2005h and spectroscopicallv ([chapman et al.ll2005h . In this 
work we intend to expand on previous results by includ- 
ing new Subaru photometry (extinction corre cted as for the 
previous dwarfs, with E(B-V)=0.08 from the lSchleeel' et al.l 
l|l998h maps), an d update the kinemat ic results obtained 
with DEIMOS in IChapman et al.l l|2005l ). having corrected 
for the error in the initial geometric models as discussed in 
§ 2. We present our observations of And IX in the same for- 
mat as the other M31 dwarfs, and list all observed stars in 
Table [5] 



3.4-1 Photometric properties 

Using the Subaru data for the satellite, the centre of And IX 
is derived as (a,5) = (00 ,l 52 m 51.1 s , +43° 11' 48.6") using the 
same technique as for the previous dSphs. In Fig. I14f a) we 
show the Subaru CMD for all stars within a 2.5' radius of the 
centre of AndIX(pink squares represent the DEIMOS can- 
didate members in this data set). The background-corrected 
radial profile was calculated using circular annuli centred on 
these coordinates and can be seen in Figure [14] (with a back- 
ground level of 23±4.8 stars per square arc minute). King 
(with r t = 6!6 ± 0.3 ), Plummer and exponential models 
have been plotted over the data. These yield half-light radii 
of 2'.7 ± 0.2, 2'.Q ± 0.1 and 2?5 ±0.1 respectively, with the 
exponential profile showing the best fit to the data (reduced 
chi-squared of 1.8 vs. 1.9 and 2.0 for King and Plummer 
profiles). Ag ain, this is con s istent with the results of the ML 
algorithm of iMartin et al.l |2008h (N. Martin, private com- 
munication) . 

Adopting the same TRGB Monte Carlo procedure as 
with And XI, we find D maa; =770 kpc for the satellite, 
and we calculate a range of distances of 615-770 kpc for 
the dwarf. From the luminosity function, we infer a dis- 
tance modulus based on the HB of 24.4 (760 kpc), also 
sh own in Figure | 14[ which is consistent with the findings 
of iMcConnachie etal.1 (|2005l ). For this reason we use their 
distance of D=765+^ 50 kpc for the remainder of our analy- 
sis. We calculate a total luminosity for the satellite at this 
distance of L = 1.491^;^ x 1O 5 L0. These revised distance 
and luminosity values update the derived parameters of 
And IX with M„ becoming -8. lip'* (previously estimated as 
-8.3), and using the result from the Plumm er profile, r 1 / 2 = 
552+ ^n pc, previo usly estimated as 530 pc (|Chapman et al.l 
l2005h and 300 pc l|Harbeck et al.ll2005h . 

As for the previous dSphs, we estimate t he mean pho- 
tomet ric metallicity of And IX by overlaying iDotter et al.l 
(2008) isochrones (with age and abundance as detailed in 



§ 3.1) for the range of distances calculated above (615-770 
kpc), giving a metallicity range for AndIX of [Fe/H]=-1.9 
to -2.3. We overlay the best observed fit for the HB distance 
(765 kpc) in Fig[Tl with [Fe/H]=-2.2. 

3.4-2 Spectroscopic properties 

We display our kinematic results for AndIX in Fig. 1151 
As with And XIII, the third line of the Ca II triplet in 
the stars observed in the AndIX masks experiences signif- 
icant contamination from OH sky lines, so we derive ve- 
locities using the first 2 lines only. We identify nine stars 
as potential members of the satellite, at a velocity of ~- 
200 kms" 1 (highlighted by the heavy histogram) and an un- 
corrected a v =5.9 kms -1 . Like And XIII, the systemic veloc- 
ity of And IX places it much closer to the regime of Galactic 
foreground contamination. AndIX also sits much closer (in 
projection) to M31 itself, so contamination from M31 halo 
stars is a larger problem than for the other satellites, due to 
the increased density in this region. Given both these fac- 
tors, it is important for us to fully check our sample for any 
contaminants. One star (highlighted in pink) sits at a dis- 
tance of roughly 10' from the centre of the satellite, which 
is beyond the tidal radius of AndIX(6.6'). Further, when 
checking its position in the Subaru CMD and its spectro- 
scopic metallicity when compared with the other eight can- 
didates, we find that this anomalous 9th star lies well off 
the RGB for the satellite (highlighted by the green circle in 
Fig. I14[) and is far more metal rich spectroscopically when 
compared to the remaining eight stars ([Fe/H]=-0.6 versus 
[Fe/H]~-2 to -3 for the other candidates). For these reasons, 
we exclude this 9th star from the remainder of our analysis. 
Two other stars with velocities consistent with the systemic 
velocity of AndIX at -220.9 kms" 1 and -209.8 kms" 1 have 
very low TDR values of 2.7 and 1.6 respectively, falling well 
below our quality cut of TDR>3.5, indicating that the ve- 
locities for these objects are unreliable. Therefore, we do not 
include these objects in the rest of our analysis, and label 
them as tentative members of And IX. 

We display the spectra for each of the remaining 8 stars 
in Fig. [15] also, with the top panel displaying the summed 
spectrum for the satellite. As was the case for And XIII, it 
can be seen that the CaIIg662 is obscured by OH skylines for 
the majority of members. This is also evident in the com- 
posite spectrum, where the equivalent width of the third 
line is measured to be significantly smaller than that of the 
first line (whereas in And XIII, the EW of the 3rd line was 
enhanced). If we were to include this line in our analysis, 
it would cause us to underestimate the value of [Fe/H] for 
the dwarf. For this reason, we derive the metallicity for the 
satellite based on the CaIIg4gg and CaIIg542 lines only, giv- 
ing a metallicity for AndIX of |Fe/Hl=-2.2 ±0.2, which is a 
more metal poor value than that found bv lChapman et all 
(2005) of [Fc/H]=-1.5, and more consi s tent w ith the value 
of [Fe/H]=-2.0 found bv lHarbeck et al.l (|2005l ) and the pho- 
tometric metallicities we quote above. We shall discuss this 
discrepancy further in § 4. Such a low metallicity would make 
this dwarf one of the most met al-poor M31 dSphs o bserved 
to date, comparable to And V (|Davidge et al.|[2002h . 

A maximum likelihood distribution for And IX is shown 
in Fig. [16] showing the best value when correcting from the 
instrumental measurements. From this technique we resolve 
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Table 5. Properties of candidate member stars in AndIX, centered at «=00:52:51.1, <5=+43:ll:48.6. The two stars with unreliable 
velocities (TDR< 3.5 and the contaminant star at >rt are also included as the last three entries. 



star a 


star 8 


vel (kms 1 ) 


[Fe/H] spec 


[Fe/H] phot 


S:N ( A" 1 ) 


TDR 


/-mag 


V-mag 


00:52:45.74 


+43:12:53.6 


-198.8+5.8 


-2.10 


-2.77 


2.6 


4.5 


20.8 


21.9 


00:52:53.18 


+43 12 48.2 


-208.3+2.4 


-2.20 


-2.75 


4.7 


7.5 


20.9 


22.1 


00:52:58.04 


+43:11:40.6 


-201.2+5.1 


-2.34 


-2.04 


2.7 


4.1 


20.9 


22.2 


00:52:58.96 


+43:13:39.1 


-215.0+4.0 


-1.90 


-2.18 


2.5 


4.3 


21.0 


22.2 


00:52:53.35 


+43:12:41.3 


-197.7+6.0 


-2.00 


-2.55 


3.1 


8.3 


21.1 


22.2 


00:52:52.30 


+43:12:41.1 


-210.7+2.3 


-1.75 


-1.73 


3.7 


12.1 


20.9 


22.3 


00:53:15.36 


+43:14:11.1 


-209.3+12.1 


-2.84 


-2.98 


3.2 


2.7 


21.5 


22.5 


00:52:50.32 


+43 11 54.4 


-201.7+20.3 


-2.90 


-3.00 


2.5 


1.6 


21.5 


22.6 


00:53:49.38 


+43 12 20.0 


-202.3+4.5 


-0.60 


-0.61 


3.1 


5.0 


21.6 


23.9 



tV^-207.7^2'2 kms -1 and a v — 4.5I3 4 kms -1 (overlaid as a 
Gaussian curve in the top panel of Fig. 1151) . T his dispersion is 
significantly lower than the value quoted in lChapman et ahl 
(2005) of a v =6.8 kms -1 , and this is likely due to the in- 
clusion of non-members in their analysis as a result of the 
incorrect geometric model. 

As for And XII, this measured velocity dispersion al- 
lows us to constrain the dark matter content of AndIX 
(rather than just assign an upper limit, as for And XI and 
And XIII), assuming virial equilibrium. As before, we first 
use eqn. (2) to calculate the M/L of AndIX. We derive 
fio =28.42 mags/arcsec 2 and So=0. 155+0. 01 Lq/pc 2 , and 
take ao = 4.5^3'! kms -1 from our ML analysis, giving a 
most probable mass to light ratio M/L = 78± 1 g°M Q /Z/e. 
The implied mass is therefore 1.2j^'| x 10 7 M©, potentially 
low mass when compared to other bri ght dwarfs (My < —8). 
Once again, applying eqn. (3) from llllingworthl (| 19761 ). we 
calculate a mass of 9.9tg 2 9 6 x 10 6 M Q for And XI, with a 
corresponding M/L=66, which again highlights the incon- 
sistencies between the two methods. 

Finally, applying equation (4) to our data, we estimate 
the mass contained within the half-light radius of And IX to 
beM hai/ = 6.5lg;5xlO 6 M . As these 6 And IX members sit 
within the half-light radius we determined in § 3.4.1, we can 
consider our velocity dispersion to be a good approximation 
for the central dispersion of the dwarf. 



4 DISCUSSION 

4.1 Comparison with previous studies 

The 4 dSphs analysed in this wor k have all been pre- 
viously examined photometrical l y ([Harbeck et al.l 120051 ; 
IZucker et all 120041 : Wartin et all 120061 ) a nd, in the case 
of AndIX and And XII, spectroscopically ijChapman et al.l 
I2005l . l2007l ). Here we discuss how our results compare with 
these other studies. For Andromeda's XI, XII and XIII, the 
deeper photometric images from Subaru, combined with our 
spectroscopic results from DEIMOS, have enabled us to sig- 
nificantly improve upon previous values for structural pa- 
ra meters, such as half-light radius, distance etc., as derived 
in lMartin et~aH (|2006h from CFHT-MegaCam data. We dis- 
play our results for the structural properties of the dwarfs 
in table [6] and their kinematic properties are displayed in 
table [7] Whilst the absolute magnitudes mirror the initial 



values very closely, we find the other structural parameters 
for the satellites all differ significantly. In all three cases 
the half-light radii have increased from ~115 pc to ~200 
pc which is comparable to what is observed for MW dwarf 
galaxies of a si milar luminosity, e .g. LeoT with My=-7 . 1 and 
r 1/2 = 170±15 (jlrwin et al.ll2007l : ISimon fc Gehall2007f ). This 
is in contrast to the brighter M31 dSphs (My <-8), which 
are all observed to be larger than their MW co unterparts 
of a similar luminosity l|McConnachie et al.ll2005h . We have 
also been able to derive distances to each of the dwarfs indi- 
vidually, whereas in Martin et al. (2006), photometry would 
only allow for an ensemble average. We find And XI to be 
the closest dwarf to us, lying at a distance of 760 kpc, with 
both And XII and And XIII being more distant, located at 
>800 kpc. In addition, we have reanalysed And XII spectro- 
sc opically, and our findin gs are largely consistent with those 
of IChapman et all (|2007l ). 



For Andromeda IX, the structural properties we have 
derived from the Sub ar u da ta closely mirr o r the find- 
ings of IHarbeck et all (|2005l ). IZucker et all l|2004h and 
IChapman et al.l (|2005h . Our range of distances (615-770 kpc) 



and our HB distanc e estim ate of 760 kpc agree well with the 
iMcConnachie et all l|2005l ) result of 765 kpc (using a TRGB 
method). Our value for the half-light radius of 552^ 22 n pc 



is sign ificantly larger than that quoted in Harbeck et al 
12005), but agrees with the findings of IChapman et al 



( 2006h and, like the remainder of its bright M31 brethren, 
it is more extended than its Milky Way counterparts of a 
si milar luminosity (e.g . Draco, with My=-8.8 and ^2=221 
pc. lMartin et al.ll2008l ). However our re-reduction and anal- 
ysis of the spectroscopic data with the correct geometrical 
model causes us to update a numbe r of th e parameters orig- 
inally presented in I Chapman et all (|2005l ). In our analysis, 
we resolve a velocity dispersion of only 4.5I34 kms -1 , com- 
pared with their measured value of 6.8 kms -1 . Reassessing 
the combined spectrum for AndIX, we find [Fe/H]=-2.2±0.2 
for the satellite, significantly mor e meta l poor than [Fe/H]=- 
1.5 reported by IChapman et al] |2005), largely due to ex- 
cluding the CaIIg662 line from our analysis, but also a result 
of the misclassification of member stars in IChapman et al.l 
l|2005l ). This new value i s in m uch better agreement wit h the 
work of IHarbeck et al.l (|2005l ) and lZucker et all (|2004l ) and 
is supported by our deep Subaru photometry. 
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Table 6. Structural properties of And XI, And XII, And XIII and And IX as derived in this work 



Property And XI And XII And XIII And IX 



a h 46 m 21 a Q h At m 2T h 51 m 51 s h 52 m 51.1 s 

<5 +33° 48'22" +34° 22'29" +33° 00'16" +43° 11'48.6' 

My,o -6-9±g;f -6.4±°;i -6.7±° : ? -S.llg;* 

r V 2 145l^ 2891™ pc 203^7 pc 5521% pc 

Distance 760±ig k P c 830 ±30° k P c 910 ±160 k P c 765 ±150 k P c 



Table 7. Kinematic properties of And XI, And XII, And XIII and And IX as derived in this work 



Property And XI And XII And XIII And IX 

v r (kms- 1 ) -419. -558.4±^j -195.0±g4 -207.7^ 2 2 7 

a v (kms" 1 4.6* " 2.6^^' 9.7±f:§' ^tf 6 4 

[Fe/H] spec -2.0±0.2 -2.1±0.2 -1.9±0.2 -2.2±0.2 

[Fe/H] phot -2.0 -1.9 -2.0 -2.2 

* Value quoted represents upper limit set by la uncertainty from maximum likelihood analysis 




Figure 17. Mean values for [Fe/H] for both M31 and MW 
satellites vs. absolute magni tude (My). A lin e of best fit to 
the MW dwarfs, taken from iKirbv et al. | J2008T) is overplotted. 
Values for the M31 dwarfs were 



12002], 


l2000t); iDavidge et al.l 


1 1999I1: 


iMcConnachie et al.l 


(2006a; 


; McConnachie et al.l 




iGrebel &: Guhathakurtal 
[M^C^rm^chie^TTroin] 
; llrwin et al. l l200g|) ; 



(2006,2007a , 2007b)TMartin et al] [12007! ); Ilbata et all &200dt) ; 
IKirbv et aD l|2008h : iMateol l|l998r i. MW objects are plotted as 
green triangle, M31 objects as blue inverted triangles and our 
results are plotted as red squares. 



4.2 Metallicities 

The metallicity for each dwarf has been reassessed, with 
all four measured to be more metal-poor than in previous 
works. There is a good agreement between Subaru photom- 
etry (using Dartmouth isochrones) and spectroscopy for all 
the dSphs. In the past, much work has been undertaken to 
determine the metallicities of faint dSphs in the MW to see 
if they follow the same pattern of decreasing metallicity with 



decr easing lumin osity as seen in their brighter counterparts 
(e.g. lMateolll998h . Until recently, this trend has appeared to 
break down in the fainter regime of dSphs, with low lumi- 
nosity satellites measured to be significantly mo re metal-rich 
than expected. However, work undertaken by IKirbv et al.l 

(2008) , measuring iron abundances with high and moder- 
ate resolution spectra from the HIRES and DEIMOS in- 
struments on Keck, has shown that the faint satellites in 
the MW do contain sizable populations of very metal-poor 
stars. This result places them in line with the previously 
established trends, as shown in Figure 1171 where we plot 
the relation between absolute magnitude, M„ and [Fe/H] 
for dSphs. The green symbols represent the MW dwarfs and 
the dashed l ine shows th e best fit to these satellites as de- 
termined by IKirbv et al.l (|2008h . We have also plotted the 
available information for the M31 dwarfs (blue inverted tri- 
angles) and it can be seen that they agree quite well with 
the relation at higher luminosities. For all four dwarfs anal- 
ysed in this paper (plotted as red squares), it can be seen 
that they appear to be consistent with the relation within 
their la error bars. This suggests that the relation between 
metallicity and absolute magnitude does not break down at 
these fainter luminosities. 

Such a finding is perhaps unexpected, given work pre- 
senting evidence for a common mass scale for the halos of 
the dSphs population. IStrigari et al] ((2008) conclude that 
the dSphs of the Local Group all have a comm on mass of 
~ lO^ Mp; within their central 300 pc. Work by I Wolf etail 

(2009) measuring the mass content of dSphs within their 
half-light radius, also presents evidence for a common mass 
scale for the halos of the dSphs population, concluding that 
even the faintest MW dSphs are consistent with having 
formed within dark matter halos of ~ 3 x 1O 9 M0. Work- 
ing within this framework of a common mass scale for the 
halos of dSphs, one might expect to see metallicity level 
off at fainter luminosities, rather than continuing to de- 
crease as, although the amount of luminous matter present 
in the system is decreasing, the amount of dark matter re- 
mains roughly constant. This would impose a limit to the 
amount of metals that could escape the gravitational poten- 
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Table 8. Mass and M/L estimates for And XI, And XII, AndXIII and And IX as derived in this work 



Property And XI And XII AndXIII And IX 



M/L flT86 (Mq/L ) 152* 

M H T86 (M s ) 7.4x10 s 
M/L /76 53* 

M J76 ( Mq) 2.6xl0 6 

M half (M Q ) 1.7xl0 6 

* Value quoted represents upper limit set by la uncertainty from maximum likelihood analysis 
RT86 refers to Richstone fc Tremaine (1986), eqn. (1) 
176 refers to Illingworth (1976), eqn. (2) 



+152 S91 +1065 117 +102 

00 -55 ozi -538 11 '-100 

1.7±^ x 10 6 3.4±^ x 10 7 1.2±J-| x 10 7 

54 + 152 401+ 505 66+ 84 

°^ -54 w -252 °°-66 

1.6±t;| x 10 6 1.6?;J x 10 7 6.5±|J x 10 6 

i'i x 106 L1 o~7 4 x lo7 * 6 - 5 ie 5 x 106 



tial of the galaxy, preventing dSphs fr om becoming overly 
metal poor. However, a similar study bv lWalker et~afl (l2009i ) 
(also measuring the mass content of dSphs within the half- 
light radius) finds that, whilst the dSphs of the MW are 
broadly consistent with having formed within a Universal 
dark matter halo, there is significant scatter in this rela- 
tionship at the faint end of the luminosity spectrum which 
could indicate that these fainter dwarf galaxies are residing 
in less massive halos. This would allow them to expel more 
of their metals at early epochs, resulting in lower metal- 
licities, which could explain this continuation of the trend 
seen in brighter dwarfs. An other explanation is provided by 
ISalvadori fc Ferraral (|2009u who model the merger history of 
the Milky Way from z=20 to the present day to investigate 
the nature of 'Ultra Faint' dSphs (such as those discussed 
in this work); they find that the low metallicities of these 
objects are the result of lower gas metallicity at the time 
of their formation, and suppression of star formation during 
their evolution. 



4.3 The dynamics of M31's dSph population 

Recent work bv lWalker et all J2009T I and I Wolf et~ai] (|2009h . 
investigating the common mass scale for dSphs demonstrates 
a clear correlation between the velocity dispersion and half- 
light radius for MW dSphs. In Fig. [18] we plot a v vs x\/% f° r 
all M W dwarfs with pub lished velocity dispersions (taken 



from Walker et al 



2009) , the availab l e M3 1 data points 



20081 ; iLetarte et all 120091 . And XVI is 



ijPefiarrubia et all 
shown as an upper limit) and o verplot the be s t fit p ower law, 
NFW and cored profiles from lWalker et aD l)2009l ) to these 
data. The upper limit determined for And XI is displayed 
as an open red square, and the resolved measurements for 
And XII, AndXIII and And IX are shown as filled squares. 
Interestingly, And XII and And IX both fall below these re- 
lations within their la uncertainties. Whilst these results 
could be an artifact of small sample size and low S:N, given 
the scarcity of available kinematic data for the dSph pop- 
ulation of M31, it is useful to consider the implications of 
such findings. If the dispersions of these dwarfs do fall signif- 
icantly below established trends, could this be an indication 
that the dSph companions of M31 are dynamically colder 
than their MW counterparts? The brighter dwarf galaxies 
of M31 are already ob served to differ from those orbiting 
the MW, with work by iMcConnachie et al.l l|2005l ) showing 
that M31 dSphs with Mv < —8 are spatially more extended 
than MW dSphs of the same luminosity (though it is un- 
clear what the origin of this difference is). Curiously, if M31 



and MW dSphs inhabit systematically similar halos (i.e. the 
mass and extent of the stellar population is independent of 
the mass and extent of the dark matter halo) , one would ex- 
pect to see the opposite trend, with the velocity dispersions 
of M31 satellites being r oughly double those of M W dwarfs 
of the same luminosity (|Penarrubia et al]|200 8). It is also 
worth menti oning that a recent study of And I, II, III, VII, 
X and XIV (|Kalirai et al.ll2009l . I2OI0I ) finds that a number 
of these objects also appear to be dynamically colder than 
their MW counterparts. 

These colder dispersions translate to lower masses for 
And IX and And XII, as the three mass estimators we have 
used within this work are all proportional to the square of 
the dispersion (as shown in eqns. (2)-(4)). Comparison of 
the masses of Local Group dSphs has previously proven to 
be problematic. A variety of mass estimators (as discussed 
in § 3.1.2) used throughout the literature, and the differ- 
ing radial extent of the dSphs probed by kinematic surveys 
make it difficult to ensure that, when making comparisons 
between dSphs within the Local Group, we are truly com- 
paring like with like. Historically, the masses of dSphs have 
been discussed in terms of the Mateo (1998) relation be- 
tween M/L and the luminosity of dwarf galaxies, which was 
established in relation to Local Group dSphs with My <-8. 
This relation is displayed in the top right panel of Fig. 1181 
with M/L values for satellites from both the MW and M31 
taken from the literature plot ted also. We display the up- 
per limits on M/L (using the iRichstone fc Tremaina Il986l 
mass estimator) for And XI (open red square), and plot the 
calculated values for And IX, And XII and And XIII (filled 
red squares), where the error bars represent their la un- 
certainties. We note that the lower limits for three of the 
four dwarfs (And XI, XII and IX) are consistent with zero. 
And IX, XI and XIII are consistent with this relation within 
their la error bars, but AndXIIjoins the ultra faint MW 
dwarfs as an outlier. 

The work bv lWalker et~ai1 (|2009l ) and lWolf et aD (|2009h 
has shown that the mass within the half-light radius of dis- 
persion supported galaxies, such as dSphs, can be accurately 
estimated with only mild assumptions about the spatial vari- 
ation of the stellar velocity dispersion anisotropy, unlike the 
two methods discussed above. Also, by considering only the 
mass contained within r 1( / 2 , one can negate the effect of the 
differing radial extents probed by surveys, making compar- 
isons between the masses of the dSphs more meaningful. 
Their conclusion that there is a tight correlation between 
velocity dispersion (and hence, mass) and half-light radius 
for MW dSphs is consistent with the idea that all the MW 
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Figure 18. For all figures, the symbols used are as for Fig 1171 In all cases, filled symbols represent resolved measureme nts, and open 
symbo ls represent upper limits set by la uncertainties. References for MW and M31 data are as specified in Fig 1171 plus lWalker et abl 
l|2009l l and references therein. Top left Global velocity dispersion for MW and M31 dSphs are plotted as a function of t heir half-light 
radius , The best-fitting mass profiles (based on power-law, NFW and cored halo profiles) for the MW dSph population, from lWalker et all 
(2009) (transferred into the a v -r 1 / 2 plane) are overplotted. Top right: Comparison of the mass-to-light ratios (M/L) and luminosity 
estimates (M„) for the faintest known dwarf galaxies from both the MW and M31 with central velocity dispersion estimates. The Mateo 
(1998) relations for dwarf spheroidals are over plotted as black lines, representing curves of constant dark matter halo mass (1, 2, 4, 8x 
10 7 Mq from bottom to top), assuming a stellar mass-to-light ratio of 2.5 Mq/Lq. And XII appears to be an outlier to this relation. 
Bottom left: Comparison of the mass contained within the half-light radius {M-half) as a function of rjyj, with the IWalker etajj 
fcoogl ) mass profiles overplotted. Both And IX and AndXII fall below these relations. Bottom right: Mean density (< p >) within the 
half-light radius is plotted, with the same mass profiles as in the top right panel. Again, And IX and AndXII appear to be outliers to 
these relations. 



dSphs (including those with My > —8) are formed with a 
"Universal" halo profile (as discussed in § 4.2). In the bottom 
left panel of Fig 1181 we plot M^aif vs. t 1 / 2 for And IX, XI, 
XII and XIII along with the other Local Group dSphs for 
which kinematic data are available, and over plot the best 
fit po wer law, NFW and cored profiles from IWalker et al.l 
(2009) to these data. In this instance, we can see that And IX 
and AndXII fall below these relations within their la errors, 
whereas And XIII conforms to the MW trends. 

In terms of its mass and dispersion, And IX in particu- 
lar seems to stand out. At My = —8.1 it is similar in terms 
of its absolute brightness to the MW dwarf s, Draco and 
Ursa Minor (My = -8.8 and -8.4 respectively, iMartin et al.l 
(2008)), and has a half-light radius that is roughly twice that 



of these dwarfs (552 pc vs. ~ 200pc), however we measure 
a dispersion of only 4.5 kms~ for And IX, c ompared with 
9.1 and 9.6 kms -1 for Draco and Ursa Minor ( Walk er et al.l 
1200112009! ). This suggests that And IX inhabits a dark mat- 
ter halo that is significantly different to those of the brighter 
MW dSphs. 

In the case of AndXII, whilst it does fall below the pro- 
files of IWalker et all |2009), it is quite similar in terms of 
both structure and kinematics to several of its MW coun- 
terparts at a similar luminosity, e.g. Hercules and Leo IV 
(with My-=-6.0 and -5.1 , rx/2=330 and 116 pc, a v =3.7 and 
3.3 kms -1 respectively. I Walker et al] 1200 9). Perhaps then, 
it is only these 'Ultra Faint' dSphs (My >-8) that reside 
in significantly different dark matter halos, with And IX be- 
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ing a curious exception. It is worth noting that the upper 
limit of <t„=4.7 kms -1 that we attach to AndXI (My=-6.9, 
r!/2 = 145 pc) suggests that this faint dwarf is also an outlier 
to these relations. 

If these satellites truly fall below these relations, the 
question becomes, are they embedded in a less massive dark 
matter halos than their MW counterparts, constituting a 
population of lower mass dwarfs whose detection has thus 
far eluded us? These lower values of M/L and M^aif are 
unlikely to be the result of tidal mass stripping, as the pro- 
cess tends to increase r ather than decrease values of M/L 
ijPenarrubia et al]|2008h . This is because the tightly bound 
central dark matter "cusp" of a dwarf is more resilient to 
tidal disruption than the stellar component of the galaxy, 
causing the object to become more dark matter dominated 
as it undergoes tidal s tripping. It should be noted that 
the mass e s timat ors of iRichstone Sz Tremaind (|l986h and 
Illing worthl (| 19761 ) are measures of the central population of 
the galaxies, where the stars can be used as tracers, there- 
fore they are only true throughout the dwarf if the mass- 
to-light ratio is constant throughout the system. It could be 
that these dwarfs reside in more massive dark matter halos, 
such as is seen for brighter galaxies, and these central low- 
mass estimates could imply lower central densities for the 
two dwarfs, making them peculiar objects. To examine this 
possibility further, we plot the mean total (i.e. baryonic plus 
dark) mass density of the dSphs within the half-light radius 
in the bottom right panel of Fig 1181 and it can be seen that 
both And IX and And XII have a lower density than would 
be expected for their size. 



5 CONCLUSION 

We have studied the four faint M31 satellites, And IX, 
AndXI, And XII and And XIII, using spectroscopic data 
from Keckll/DEIMOS and Subaru/Suprime-Cam. Using 
our spectroscopic results, we have identified probable mem- 
bers for each galaxy and constrained their systemic veloci- 
ties. Owing to our low S:N and small sample sizes, we are 
unable to resolve a velocity dispersion for AndXI. Instead 
we use the la confidence level from the maximum likelihood 
technique as an upper limit for the satellite. In the case of 
And IX, AndXII and And XIII we resolve a v =4.5tl lkms~ 1 , 
o-^=3.8l3 gkms _1 and a v —9.7^t, kms -1 for each dwarf re- 
spectively, though we note that AndXII is consistent with 

kms -1 within its la errors and the measurement of 
AndXIII's dispersion is based on only three stars. Measuring 
the average spectroscopic metallicity of each satellite from 
their composite spectra, we find that all four satellites are 
metal-poor ([Fe/H]~-2.0). Such metallicities are comparable 
to what is observed in MW dSphs of a similar luminosity, 
suggesting that the established relation between My and 
[Fe/H] obser ved in both the brighter cla ssical MW dwarfs 
l|Mateoll 19981 ) and the fainter MW dSphs l|Kirbv et al.ll2008l ) 
continues down to luminosities of Mv > —6.4 in the M31 
dSphs population. 

Finally, we place the velocity dispersions and masses 
for these dSphs in the context of the universal dark mat- 
ter halo profiles est ablished for the MW dwarf galaxies in 

1 Walker et all (|2009h and find that both And IX and AndXII 
have lower dispersions than would be expected for dSphs 



of their size , simi lar to a number of dSphs studied by 
iKalirai et ail (|2010h - If these results are verified, it could im- 
ply that the dSph population of M31 is dynamically colder 
than that of the MW, which is unexpected, given their larger 
spatial extent. Such results also imply that the mass con- 
tained within the half-light radii of these dwarfs is also lower 
than what is typically observed in their MW counterparts 
of a similar size, and suggests that these dSphs may reside 
in lower density dark matter halos. 

The results we have presented here indicate that, in 
terms of their dynamical properties, these satellites at the 
faint end of the dwarf spheroidal spectrum do not behave as 
expected from previous studies. Further spectroscopic obser- 
vations of faint objects si milar to these (e.g., those found by 
iMcConnachie et al.| [2008) will undoubtedly help us to better 
understand this low luminosity regime, and shed light onto 
the behaviour of these faintest of galaxies. 
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